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Abstract 

We discuss the proton lifetime in pure gravity mediation models with non- 
universal Higgs soft masses. Pure gravity mediation offers a simple framework for 
studying SU(5) grand unified theories with a split supersymmetry like spectra. We 
find that for much of the parameter space gauge coupling unification is quite good 
leading to rather long lifetimes for the proton. However, for 7713/2 ~ 60 TeV and 
tan/3 ~ 4, for which gauge coupling unification is also good, the proton lifetime is 
short enough that it could be in reach of future experiments. 


1 Introduction 


After the initial run of the LHC, the constraints on new physics are rather severe [ 1 ]. 
Although models can still be made to realize weak-scale mass spectra, sfermion masses of 
generic models like the constrained minimal supersymmetric standard model (CMSSM) 
[ 2 ] are now required to be larger than about a TeV. As a result, the naturalness of 
supersymmetry (SUSY) has been called into question. However, it was perhaps naive 
to expect nature to fall into our strict definition of naturalness with less than 10 % fine- 
tuning. Supersymmetry, with sfermion masses larger than a TeV, still solves the larger 
hierarchy problem associated with grand unification and/or the Plank scale. Furthermore, 
if the sfermion masses are set by the gravitino mass, 777.3/2, and are larger than about 10 
TeV, the gravitino lifetime is short enough that it decays before BBN [ 3 ]. Moreover, as 
the mass scale of the sfermions is pushed beyond the weak scale, the constraints on SUSY 
models from flavor and CP violation in the sfermion sector are greatly relaxed [ 4 ]. These 
advantages, plus the fact that sfermion masses this large are consistent with a larger 
Higgs mass like the 126 GeV Higgs boson seen at the LHC [ 5 ] suggest we relax our strict 
definition of naturalness. 

Large sfermion masses like those found in split supersymmetry [ 6 ] are realized in 
models such as pure gravity mediation (PGM) [ 7 , 8 ], which can be parametrized by a single 
parameter [ 9 ] 7773/2- This minimal model of pure gravity mediation is similar in many ways 
to minimal supergravity (mSUGRA). Universal masses equal to 7773/2 are imposed at the 
grand unified theory (GUT) scale based on the assumption that the Kahler manifold is 
flat for all matter fields. Unlike the CMSSM, gauginos do not get a tree-level mass. This is 
because the supersymmetry breaking field is not a singlet and so is excluded from coupling 
to the gauge kinetic function to leading order. Thus, the leading order contribution to the 
mass of the gauginos comes from anomaly mediation [ 10 ] and is loop suppressed relative to 
the sfermion masses. The U-terrn, which contributes to electroweak symmetry breaking, 
is identical to that in mSUGRA, B = A — 7773/2- However, since the A-terms of PGM are 
effectively zero, B = —7773/2 and B is fixed for a given value of 7773/2. This makes radiative 
electroweak symmetry breaking (EWSB) difficult. However, by adding a Giudice-Masiero 
term [ 11 ], B is no longer fixed by 7773/2 alone, but also depends on the coupling of the 
Giudice-Masiero term. This additional freedom in B makes radiative EWSB possible, but 
only for small values of tan / 3 . Once the Higgs mass constraint is taken into consideration, 
these models have a single free parameter which is some combination of 7773/2 and tan /3 
[ 9 ]. However, because tan /3 is restricted to be less than about 3 , 7773/2 tends to be rather 
large. The constraints on tan /3 can be removed, if the Higgs soft masses at the GUT scale 
are taken to be non-universal [ 12 ]. In this case, 7773/2 can be taken to be smaller for larger 
values of tan / 3 . 

Another important motivation for SUSY is grand unification [ 13 ]. In the Standard 
Model (SM), the gauge couplings approach each other as they are run up to the high scale 
[ 14 ]. However, the quality of the coupling unification is less than convincing. If the SM is 
supersymmetrized, on the other hand, the unification of the gauge couplings becomes quite 
good [ 15 ]. Furthermore, grand unification in the SM would generate enormous quadratic 


1 


divergences for the Higgs boson. However, these quadratic divergences are significantly 
reduced for supersymmetric grand unified theories, even if the sfermions are larger than 
a TeV. Clearly, grand unification is another motivation for PGM. 

The signatures of these simple PGM-type models are limited. One possible signature 
at the LHC for small 777.3/2 is the wino [16]. For larger 7773/2, on the other hand, the wino 
cannot be seen at the LHC but could be a viable thermal relic dark matter candidate 
[17]. If this is indeed the case, it could be seen by indirect detection experiments in the 
near future [18]. However, this scenario is already under tension from existing indirect 
detection experiments [19]. The direct detection of wino dark matter is challenging as its 
scattering cross section with a nucleon is as small as 10 -47 cm 2 [20]. A Higgsino signature 
at the LHC is another possible observable which arises from tuning fi to be small [21]. 
However, this is also difficult to see. This scenario could also have Higgsino-like dark 
matter which could possibly be seen in future indirect detection experiments [22], The 
scattering cross section of the Higgsino with a nucleon is dependent on the size of the 
wino component of the LSP, and may be probed in future experiments [23]. 

In this work, we will examine another possible signature of these models. Since the 
colored triplet Higgs gives threshold corrections to the gauge couplings when integrated 
out, the quality of the coupling unification determines the mass of the colored triplet 
Higgs [24, 25, 26] and so affects the lifetime of the proton. When the colored Higgs 
is integrated out, it also generates a dimension-five operator proportional to down-type 
Yukawa couplings which lead to proton decay [27]. Since this dimension-five operator 
is proportional to the down-type Yukawa couplings, it will be enhanced for large tan /3. 
Proton decay from this dimension-five operator arises from a loop diagram with a Higgsino 
mass insertion [28]. Proton decay of this type can then be suppressed for small /j. When 
unification is not ideal and tan /? is large, a larger Higgsino mass can increase the rate of 
proton decay from this dimension 5 operator. Parameters of this size are viable in PGM 
models. Since proton decay of this type is also suppressed by 7773/2, the more interesting 
parameter space will be for smaller 7773/2 and larger tan/L Therefore, we will need to 
consider non-universal Higgs soft masses. We will find that if 7773/2 is small and tan/3 
is larger, which is also consistent with the Higgs mass measurement, the proton lifetime 
may be in reach of future experiments. However, for much of the parameter space the 
lifetime tends to be well beyond the reach of future experiments. We will also look at 
the quality of the gauge coupling unification determined by the deviation of the colored 
triplet Higgs mass, Mh c , from the GUT-scale as well as the deviation of (ilL^Ms) 1 ^ 3 , 
where X represents the GUT scale SU(5) gauge bosons that become massive and £ is the 
24 which breaks SU(5) at the GUT scale. 

2 Mi nim al SUSY SU(5) GUT 

In this section, we will outline the SU(5) SUSY GUT theory [29, 30] we will consider. 
Additional details on these models can be found in Appendix A. The superpotential for 
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this minimal SU(5) SUSY GUT is given by 


II — Wffiggs T fUVukawa j (1) 

where 

lUffiggs = ^AyTrE 3 + ^m s TrE 2 + \ h HT,H + m H HH , (2) 

o z 

WVk.™ = - x/2/«<lf ^Hi , (3) 

and a,b , ••• = 1-5 represent the SU(5) indices and is the totally antisymmetric 

tensor with 612345 = 1. dh and T,; are the chiral superfields in the 5 and 10 representations, 
respectively, with i denoting the generation index. H and H are the 5 and 5 containing 
the mininal supersymmetric Standard Model (MSSM) doublets. In these expressions, we 
have assumed i?-parity conservation which forbids terms like TUM* and H§. The adjoint 
Higgs field, E, gets a vacuum expectation value (VEV) in the direction 

(E) = V ■ diag(2, 2, 2, —3, —3) , (4) 

breaking the SU(5) gauge group to the SM gauge groups SU(3)c'<8>SU(2)£®U(l)y. Be¬ 
cause SUSY remains unbroken for SU(5) breaking, we have V = my/Ay. For this setup, 
the masses of S 3 , E 8 , E 2 4 , and H c are given as 

My = My g = My 3 = ^AyU , My 24 = ^AyU , M Hc = 5\ H V , (5) 

while the /i term for the MSSM Higgs fields is 

g 0 = m H - 3X h V . ( 6 ) 

As is usually done, we tune the parameter mu to realize fi o <C which is typically 
referred to as the doublet-triplet splitting . 1 In addition, the gauge interactions of the 
adjoint Higgs field yield an Y-boson mass of Mx = hV2g 5 V where g$ is the unified gauge 
coupling constant. The components E( 3 * 2 ) and S ( 3 2 ) become the longitudinal component 
of the X bosons, and thus do not appear as physical states. 

The Yukawa couplings and f lJ in Eq. (3) have redundant degrees of freedom, most 
of which are eliminated by the field redefinition of T and $. Since h lJ is a symmetric 
matrix, hU and have six and nine complex degrees of freedom, respectively. The field 
redefinition of the SM fields forms the U(3)®U(3) transformation group, and thus the 
physical degrees of freedom turn out to be (12 + 18) — 9 x 2 = 12. Among these degrees 

1 This is another fine-tuning besides that for the Higgs mass. Note that the fine-tuning for the Higgs 

mass becomes worse as the g parameter is taken to be larger, while the doublet-triplet fine-tuning becomes 

less severe. This tension may explain why g is much larger than the clectroweak scale [31, 32]. 
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of freedom, six of them are the quark mass eigenvalues and four are for the Cabibbo- 
Kobayashi-Maskawa (CKM) matrix elements and we arc left with two phases [33]. In this 
paper, we take the same basis used in Ref. [25] such that 

h ij = e^Sijfu^Qc) , ( 7 ) 

f ij = V^f dj (Q G ) , ( 8 ) 

where f Ui (QG) and /<f. ( Qg ) are the up-type and down-type Yukawa couplings, respectively, 
at a scale Qg around the GUT scale, and V t j is the CKM matrix. The phase factors ipi 
satisfy the condition JT (p^ = 0, and thus only two of them are independent. In this basis, 
the MSSM superfields are embedded into the SU(5) matter multiplets as 

dh 9 {Qi, e-^Ui, VijEj} , d>i 9 (D h Li} . (9) 

Then, Eq. (3) leads to 

^Yukawa = fuAQi ' H 2 )U ia - V*f dj (Q° ■ HQD ja - f • HQ 

- ■ Qi)H‘c + V-jfPQ? ■ Lj)~Hc a 

+ - V-f^e-^e^UivDjJJcc . ( 10 ) 

The new phase factors appear only in the couplings of the color-triplet Higgs multiplets. 


3 Mass Spectrum and Coupling Unification 


To compute the proton decay rate, we need to evaluate the masses of the GUT-scale 
particles which induce the baryon-number violating interactions. In this section, we es¬ 
timate these masses using the method discussed in Refs. [24, 25, 26]. The mass of the 
heavy particles is determined by first RG running the couplings to the scale where they 
approximately unify. Then, because the thresholds at the GUT scale depend on these 
superheavy particles, their masses can be determined by assuming the deviation in gauge 
coupling unification is solely due to these thresholds. Note, we will use the DR scheme [34] 
in the following calculation. At the scale Qg near the GUT scale, the one-loop matching 
conditions for the gauge coupling constants are as follows [35, 36]: 


9 2 AQg ) 
1 
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where g G is the unified gauge coupling constant. Note that the conditions do not include 
constant (scale independent) terms since we use the DR scheme for renormalization. 
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Assuming the above equations contain the major thresholds for the gauge couplings, they 
can be used to solve for the masses 


3 2 1 _ 3 / Q g \ 

£(Qg) 93(Qg) 9i(Qg) IOtt 2 U \M Hc J ’ 

_5_3_2_ = __9_ (9± 

9i(Qg) 9i (Qg) di (Qg) 27t 2 \M g 

with Mq = The above expressions allow us to find the masses of the heavy 

particles in the combination 2 , Mh c and M 2 -My;. The value of Mh c and Mq found from 
these relationships will be used below to find the lifetime of the proton. 

4 Proton Decay 

In the minimal SUSY GUT, proton decay is induced by the exchange of the color-triplet 
Higgs boson, and the dominant decay mode is, generally, p —> K + v [27]. We will only 
give details of the contributions from the colored Higgs boson since it will often be the 
dominant source of proton decay in PGM. At the GUT scale, the triplet Higgs boson 
is integrated out. The most important interaction for our considerations is the color- 
triplet. Higgs exchange which we match at the scale Qg on to the dimension-five effective 
Lagrangia.n 

cf = + c«“o“ + h.c., ( 13 ) 

where the effective operators O r -- J kl and O r -Jh are defined by 

o% = J d 2 e i ■ Q-XQi ■ U ), 

oy, = f (Po e ahc UiXEjUkbDi c , (14) 

and the Wilson coefficients Cl^f and C^ 1 are given by 

df (Qa) = , 

M Hc 

C$?(Qa) = T )~A,V V V; ! Ue-‘»‘ . (15) 

Note, the color indices must be completely antisymmetric for these interactions and as a 
result, only operators with at least two generations will be allowed. For this reason, the 
dominant decay modes contain a strange quark in their final state, he., p K + v. 

As can be seen in Eq. (10), at the GUT scale the lepton and down-type quark Yukawa 
couplings should be equal. However, in running up from the weak scale, we find them to 

2 The third condition is used to determine ^(Qg)- 
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be quite different especially those for the first two generations. The difference is, however, 
easily compensated by effects above the GUT scale; for instance, the higher-dimensional 
operators induced at the Planck scale contribute to the Yukawa couplings, which may 
account for this difference [37, 38, 39]. Because it is not known which of these values is 
close to the correct value for the Yukawa coupling at the GUT scale, in the the discussion 
below, we use both the down quark and lepton-type Yukawa couplings to calculate the 
proton lifetime. This will allow us to quantify our uncertainty in the lifetime of the proton. 

The relevant operators in Eq. (15) can be further reduced by keeping only those with 
the largest Yukawa couplings. We find that only the operators 0 ^ 12 an d O 3311 yield a 
sizable contribution to proton decay, even though the contribution is suppressed by a 
flavor changing element of the CKM matrix. This contribution turns out to be dominant 
because of the large third generation Yukawa couplings involved [28]. The relevant Wilson 
coefficients are then 

C 5 T«g) = -j^-f,UQa)e- in V ti 1C, , 
m h c 

cIrHQg) = tt -fMQc)e-‘*'v a v;,. (ie) 

Mhc 

Notice that the coefficients include a common phase factor which is therefore not 

important for proton decay. 

The Wilson coefficients in Eq. (16) are then evolved down to the SUSY scale. At the 
SUSY scale, the sfermions of these dimension-five operators are integrated out via the 
one-loop diagram found in Fig. 7 of Appendix B. The process proceeds via the exchange 
of either a charged wino or a Higgsino. 3 In PGM, we generally have |/i| [il^l and so 
the contribution from Higgsino exchange dominates [41]. 4 For these reasons, we focus on 
the charged Higgsino exchange process in what follows. 

The loop diagram in Fig. 7 is then matched onto the baryon-number violating four- 
fermion operators [42, 43, 44] 

= c i e abc (u a m d b Ri )(Q c L3 ■ L L3 ) , (17) 

with 

C,(Qs) = ^CS U (Qs)F(^mf it ,m 2 rit ) , (18) 

where i = 1,2, and Qs is the SUSY breaking scale taken to be around m 3 / 2 . The loop 
function F is found in Appendix B. The above expression shows that the proton decay 
rate depends on the SUSY spectra through the loop function. We will see this dependence 

3 This is the dominant contribution to proton decay, unless there is flavor violation in the sfermion 
sector. In this paper, we assume there is no flavor violation in the sfermion sector. The flavor violating 
case is discussed in Ref. [40]. 

4 Higgsino exchange dominates in this limit because the gauginos and Higgsinos in the one-loop 
diagrams are required to flip their chirality, and thus their contribution to proton decay is proportional 
to their masses, as can be seen from the expression for the function F given in Eq. (41). 
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in Sec. 6 for the PGM scenario. Note that the loop function is suppressed by the sfermion 
masses. Thus, we expect that for large m 3 / 2 the proton lifetime is long enough [41, 45] 
to evade the current bound, r(p —> K + V ) > 5.9 x 10 33 years [46]. This can be compared 
to the weak-scale SUSY scenarios; in these cases, the proton decay rate is in general 
predicted to be so large that the minimal SUSY GUT is excluded [47] and thus some 
additional conspiracy is required to realize a SUSY GUT. 

We now run the Wilson coefficients down to the hadronic scale, Qhad = 2 GeV. The 
Lagrangian at this scale takes the form 5 

£{p -A K + V T ) = C usd [e abc (u a R s b R )(d c L v T )} + C uds [e abc (u a R d b R )(s c L u T )\ . (19) 

Using these Wilson coefficients, we then evaluate the partial decay width of the p —> K + V 
and fold 

T(P K+V) = (1 - 2 \ A{p _> K+V) | 2 , (20) 

327T \ m* J 

where m p and rriK are the masses of proton and kaon, respectively, and 

A(p -A K + V ) = C usd (Q had )(K + \(us) R d L \p) + C uds (Q had )(K + \(ud) R s L \p) . (21) 

The hadron matrix elements in the above equation have been recently computed in 
Ref. [48] using a lattice simulation of QCD, 

(K + \(us) R d L \p) = -0.054(H)(9) GeV 2 , 

(K+\(ud) R s L \p) = -0.093(24)(18) GeV 2 , (22) 

where the first and second parentheses represent statistical and systematic errors, respec¬ 
tively. The matrix elements are computed at the scale Qha,d = 2 GeV. 

Before concluding this section, we comment on other possible contributions to proton 
decay. Firstly, the dimension-five baryon-number violating operators in Eq. (14) can also 
be generated at the Planck scale, Mp. If the coefficients of the operators are 0(1/Mp), 
that is, there is no suppression from Yukawa couplings, then they will give the dominant 
contribution to proton decay and result in a lifetime which is too short [49]. It is expected, 
however, that there is some underlying mechanism such as a flavor symmetry which is re¬ 
sponsible for the structure of the Yukawa couplings. This symmetry could give additional 
suppression to these Planck-scale operators. In this paper, we assume that the contribu¬ 
tion of these operators is less significant compared with the color Higgs contribution, and 
neglect them in the following analysis. 

Secondly, the exchange of the X bosons will also induce proton decay. This decay 
mode is via a dimension-six GUT-scale effective operator and is thus usually subdomi¬ 
nant compared to the contribution of the dimension-five operator discussed above. An 
approximate expression for the lifetime of the proton from the dimension-six operator is 

rHaV),3xl0^x( 1Cx * GeV ) 4 . (23) 

5 For more details of how we arrived at this expression see Appendix B. 
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There is a slight dependence on the masses of SUSY particles we have neglected. As can 
be seen from this expression, the proton decay width from the dimension-six operator will 
in general give lifetimes too long to be detected, at least much longer than the present 
bound: r(p —* e + 7r°) > 1.4 x 10 34 years [50, 51]. 

5 Pure Gravity Mediation 

As discussed above, the lifetime of the proton depends on the SUSY parameters. Mo¬ 
tivated by the 126 GeV Higgs boson [5] and other cosmological considerations [52], we 
will analyze the proton lifetime for PGM models. The scalar potential of PGM takes the 
same form as that of mSUGRA 

2 

+ (A 0 W (3) + H 0 fU (2) + h.c.) + , (24) 

which is determined by the flat Kahler manifold 6 and the superpotential W is given in 
Eq. (1). HU 2 ) and are the bi- and trilinear parts of the superpotential. For PGM, 
the SUSY breaking field is a non-singlet and strongly stabilized [54] which suppresses 
the gaugino masses and A-terms respectively. The gaugino masses are regenerated by 
anomalies and take the form 7 


„, 33 g\ 

Ml = 5 16^“ 3/2 ’ 

(25) 

^ ~ 16^ ms/2 ’ 

(26) 

g\ 

M3 = 3 167T 2m3/2 ' 

(27) 

In order to account for radiative EWSB, mSUGRA is further modified by including a 

Giudice-Masiero term for the Higgs fields in the Kahler manifold [11]. 
Higgs boson parameters to 

This modifies the 

P = PO + CH m 3/2 j 

(28) 

Up = po(A 0 - m 3/2 ) + 2c H ml /2 , 

(29) 


V = 


dW 


dcj) 1 


where po is the superpotential Higgs bilinear term found in W^ 2 \ This allows us to vary 
both p, and Up independently in order to satisfy the EWSB conditions. This leaves 777 . 3/2 
and ch as free parameters. In this case, tan j3 is an output of the EWSB conditions, but 
in practice one can trade ch for tan/3 and use m 3 / 2 and tan fj as free inputs. Since this 
simplest of PGM models tends to require small tan/3 and larger 7713 / 2 , we will allow the 

6 If the Kahler manifold for the first two generations is no-scale like, these models can explain g — 2 
experiments [53]. However, in this case the proton decay calculation is more complicated because of an 
additional wino contribution but should give a similar order of magnitude for the proton lifetime. 

7 The A-terms are also regenerated by anomalies. However, they are too small to be of importance. 









Higgs soft masses to be free parameters. This will allow tan /3 to be larger and so allow 
for m 3/2 to be smaller [ 12 ], As was seen in the previous sections, both larger tan/3 and 
smaller m 3/ / 2 will lead to shorter lifetimes of the proton. We will not discuss the origin of 
these non-universal Higgs soft masses here. However, discussion about this can be found 
in Ref. [12], Lastly, we note that the non-universal Higgs soft masses, mi and m 2 , can 
also be parametrized in terms of the low scale values of /i and which are otherwise 
also outputs of the EWSB conditions. We will take advantage of this in the results below 
in order to zoom in on some features of the proton lifetime. 

6 Results 

We are now in a position to discuss the proton lifetime and mass scales associated with 
gauge coupling unification in a variety of models which have varying degrees of non- 
universality in the Higgs sector. We begin by displaying in Fig. 1 the mi = m 2 vs. tan/3 
plane for fixed gravitino mass. This is a one-parameter extension of the minimal (two- 
parameter) PGM model and resembles NUHM1 models [55]. In the left panel we have 
fixed m 3 / 2 = 60 TeV. For this value of the gravitino mass, the Higgs mass lies between 124 
and 128 GeV 8 for tan/3 roughly between 4-9 as shown by the red dot-dashed curves. The 
thin blue lines show the values of the LSP (wino) mass 9 and are solid for /i > 0 and dashed 
for /j, < 0. The anomaly mediated contribution to m x for m 3/ / 2 = 60 TeV is 170 GeV. At 
low tan /3, threshold corrections from the heavy Higgs bosons and the Higgsinos increase 
the mass for /i < 0 and decrease the mass for /i > 0 . At large tan/3, the wino mass, for 
both positive and negative /r, tends to its anomaly mediated value. The curves end at 
high and low values of the Higgs soft masses due to the absence of radiative electroweak 
symmetry breaking. For large and negative values of m\ = (the sign on the axis refers 
to the sign of the mass squared), the Higgs pseudoscalar mass squared is negative, and 
for large positive m\ = m 2 , the electroweak conditions yield |/i | 2 < 0 . 

The thicker black curves in Fig. 1 show the values of the proton lifetime. As discussed 
earlier, as there is some uncertainty as to how we match the Yukawa couplings at the 
GUT scale, we have results based on quark Yukawa couplings (shown by the solid curves) 
and results based on lepton Yukawa couplings (shown by the dashed curves). As one 
can see from the figure, the calculated proton lifetime is sensitive to tan /3 yet relatively 
insensitive to the value of mi ;2 for fixed gravitino mass. In general, the proton lifetime 
is lower at high tan /3 due to the increase in the down-like Yukawa couplings when tan /3 
is increased, whereas the Higgs mass increases with tan /3. For these relatively low values 
of the gravitino mass used in the left panel, the proton lifetimes based on quark Yukawas 
drop below 5 x 10 34 years only when tan /3 > 7 where > 127 GeV. The lifetime increases 
rapidly at lower tan/3 and exceeds 5 x 10 35 years when tan/3 < 4 where m/, < 124 GeV. 

8 We refer to this extended range of Higgs masses to account for the uncertainty in the calculation of 
the Higgs mass. Note also that the Higgs masses calculated here differ slightly from those calculated in 
[12] as here we are not imposing strict gauge coupling unification at the GUT scale. 

°Thc present lower limit on the wino mass from the LHC experiment is about 270 GeV [56]. 
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Figure 1: The tan/Smni^ plane for a) m 3 / 2 = 60 TeV and b) 1713/2 = 100 TeV. The 
Higgs mass is shown by the nearly horizontal thin red contours in 1 GeV intervals. The 
wino/chargino mass is shown by the thin solid (p > 0) and dashed (p < 0) contours. 
The thick black contours show the value of the proton lifetime based on the quark Yukawa 
couplings (solid) and lepton Yukawa couplings (dashed) in units of 10 35 years. Lifetime 
contours for the solid curves are labeled to the left of the contours whereas dashed contours 
are labeled to the right. 


However, the wino mass requires p < 0 and tan/3 > 6. Recall that these lifetimes are 
computed from Eq. (20) and when the lifetime exceeds 3 xlO 35 years, the dominant 
contribution to the decay rate comes from the dimension-six operator given in Eq. (23). 
Proton lifetimes based on lepton Yukawas are significantly smaller (by a factor of roughly 
20), so that Tp < 5 x 10 33 years when tan/3 > 6 and is still smaller than 2 x 10 34 years 
when tan /3 > 4. 

In the right panel of Fig. 1, we have taken 7773/2 = 100 TeV and as expected the 
Higgs mass for a given value of tan /3 is higher. The range 124 - 128 GeV now requires 
tan/3 ~ 3.5 - 6. The uncorrected wino mass is now about 290 GeV and in the figure we 
see lower (higher) wino masses when p > (<)0. The proton lifetimes are now significantly 
higher. At tan/3 = 6, the quark based value of t p determined by the dimension five 
operator is now 10 36 years and increases as tan/3 is lower. The lepton based lifetimes 
remain a factor of about 20 lower and may still be as low as 5 x 10 34 years at tan (5 = 6. 

To see more clearly the dependence of the proton lifetime on the PGM parameters, we 
show in the left panel of Fig. 2 the behavior of the proton lifetime as a function of tan/3 
for fixed m 3 / 2 = 60 and 200 TeV with mi = m 2 = 0. The lifetime falls off monotonically 
with tan /3. The ratio between the quark and lepton evaluation of r p is seen to be nearly 
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Figure 2: The dependence of the proton lifetime on tan (3 (left) for fixed m 3 / 2 = 60 and 
200 TeV, and on m 3 / 2 (right) for fixed tan/3 = 2 and 5. Both of the Higgs soft masses 
have been fixed mi = m 2 = 0. The dependence of the Higgs mass is also shown with 
its value given on the right side of each panel. The horizontal line indicates the present 
experimental bound [46]. 


constant as tan (5 is varied with a ratio of about 20. We also see the substantial increase 
in t p when m 3 / 2 is increased to 200 TeV. The experimental limit on the proton lifetime, 
r(p —> K + u) > 5.9 x 10 33 years [46], is shown by the horizontal line. In the right panel, 
we show this increase in r p with m 3 / 2 for fixed values of tan/3 = 2 and 5. In both panels, 
we also show the Higgs mass as a function of tan/3 and m 3 / 2 with its value displayed on 
the right edge of each panel. Restricting the Higgs mass to the range 124-128 GeV allows 
one to focus on the relevant ranges of either tan/3 or m 3 / 2 and hence on the predicted 
proton lifetime. 

In contrast to the proton lifetime, the relevant GUT mass scales, M# c and Mq , are 
relatively insensitive to the PGM parameter choices as seen in Fig. 3. As one can see in the 
left panel, there is very little dependence on tan/3. The mass parameter Mq = (M|Ms)s 
is always close to 10 16 GeV independent of m 3 / 2 (as also seen in the right panel). While the 
color-triplet mass is insensitive to tan /3, it does have a mild dependence on the gravitino 
mass and ranges from a few xlO 16 - few xlO 1 ' GeV. Notice that in the weak-scale SUSY 
scenario the mass of the color-triplet Higgs multiplet is predicted to be around 10 15 GeV 
[47]. A heavier color triplet mass makes the proton lifetime long enough to evade the 
current experimental bound. Furthermore, in some of the parameter space of PGM, the 
GUT-scale parameters Mh c and Mq are both of (9(10 16 ). In these cases, the threshold 
corrections at the GUT scale become very small, which implies the unification of the gauge 
couplings is quite good. In fact, for m 3 / 2 ~ 60 TeV and tan/3 ~ 5, we get good gauge 
coupling unification and a proton lifetime which could be in reach of future experiments. 

In Fig. 4, we offer two additional planes which show the dependence of the proton 
lifetime on other PGM parameters. In the left panel, we plot the lifetime contours in the 
mi = m 2 , m 3 / 2 plane. This is again a NUHMl-like model and we have fixed tan/3 = 5. 
As in Fig. 1, the red-dashed curves show the Higgs mass contours which vary from about 
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Figure 3: The dependence of the GUT-scale mass parameters (M Hc and M G ) on tan/3 
(left) for fixed m 3 / 2 = 60 and 200 TeV, and on m 3 / 2 (right) for fixed tan/3 = 2 and 5. 
Both of the Higgs soft masses have been fixed mi = m 2 = 0. The dependence of the Higgs 
mass is also shown with its value given on the right side of each panel. 


124-128 GeV for the plane shown. As before, the curves extend across a limited range in 
mi = m 2 where the EWSB conditions can be satisfied. At large positive rn\ , /i 2 goes to 0 
(where the curve is cutoff). At very small p, the Higgs masses increases rapidly causing 
the sudden downturn in the mass contours. As expected, we see the wino mass varies 
considerably as m 3 / 2 is varied. For the range in m 3 / 2 shown, the proton lifetime varies 
from as low as 10 33 years using the lepton Yukawas and low m 3 / 2 to as high as 10 37 years 
using quark Yukawas and m 3 / 2 ~ 150 TeV. 

In the right panel of Fig. 4, we show a two-parameter extension of the two-parameter 
PGM similar to the NUHM2 [57]. Results are displayed in the p,mA plane for fixed 
tan/3 = 5 and m 3 / 2 = 60 TeV. In this case, the EWSB conditions, are used to solve for 
the two Higgs soft masses which now differ. As the Higgs mass is largely independent of 
mA, the Higgs mass contours are nearly vertical. At the center of the plot, as \p\ gets to 
be very small, gets large and exceeds 130 GeV. At large \p\, is always larger than 
125 GeV in the ranges shown. The threshold corrections to the wino mass are sensitive 
to p and mA and that accounts for the variation of rn x as these parameters are varied. 

The proton lifetime varies between 10 34 and 10 36 years but shows significantly more 
variability. This is due to the competing effects of changing p. The proton lifetime 
depends both on the color-triplet Higgs mass and on p itself. 10 As p is lowered, the 
color-triplet Higgs mass decreases which tends to decrease the proton lifetime. But as \p\ 
is further decreased, the proton lifetime dependence on p overcomes its dependence on 
Mh c and the lifetime increases very rapidly at small \p\ seen by the sharp downturn in 
the contours near p = 0. These effects can be better understood by examining Figs. 5 
and 6 which show the behavior of the proton lifetime and GUT-scale masses, including 

10 The proton decay rate directly depends on p through the loop function F in Eq. (18). When 
\p\ « ?ri 3/2 , F oc /j/mU, while if p >> m 3 / 2 , F oc \og(p 2 /my 2 )/p, as can be seen from Eq. (41). 
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Figure 4: a) T/ie m 3 / 2 ^ 1,2 plane for tan /3 = 5 and b) the p^riA plane for tan /3 = 5 and 
777 . 3/2 = 60 TeV. In both panels, the Higgs mass is shown by thin red dot-dashed contours in 
1 GeV intervals. The wino/chargino mass is shown by the thin solid (p > 0) and dashed 
(u < 0) contours. The thick black contours show the value of the proton lifetime based 
on the quark Yukawa couplings (solid) and lepton Yukawa couplings (dashed) in units of 
10 35 years. Lifetime contours for the solid curves are labeled to the left of the contours 
whereas dashed contours are labeled to the right. 

the heavy Higgs mass, as a function of p for fixed tan/3 and 7773 / 2 - Here we see the first 
gradual and then rapid decrease in the color-triplet mass as \p\ is lowered from large 
values toward p = 0. There is no substantial difference in this behavior between the two 
values of tan /3 shown. Once again, we see that Mq depends very little on our parameter 
choices and is always near 10 16 GeV. 

Finally, in Fig. 6, we see the sharp increase in the proton lifetime as \p\ gets small. 11 
Here we see also that the Higgs mass rises sharply as p tends to zero. It is important 
to recall that the lifetime plotted corresponds only to that given by the dimension-five 
operator given in Eq. (20) and would not exceed 3 x 10 35 years when the dimension-six 
operator is included. The latter is fairly insensitive to parameter choices. 

7 Conclusion and Discussion 

As we await new results for physics beyond the standard model from the LHC, we have 
been forced to consider supersymmetric models with sfermion masses larger than what 

n Our calculations are only valid for \p\ much greater than the wino mass. 
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Figure 5: The dependence of the GUT-scale masses (mH c and M G ) on p for fixed 

tan/3 = 3 (left) and tan/3 = 5 (right) and fixed m 3 / 2 = 65,130, and 260 TeV. Both of the 
Higgs soft masses have been fixed mi = m 2 = 0. The extent of the curves is determined 
by the validity of radiative EWSB. 
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Figure 6: The dependence of the proton lifetime on p for fixed tan/3 = 3 (left) and 
tan/3 = 5 (right) and fixed m 3 / 2 = 65,130, and 260 TeV. Both of the Higgs soft masses 
have been fixed m\ = m 2 = 0. The extent of the curves is determined by the validity of 
radiative EWSB. The dependence of the Higgs mass is also shown with its value given on 
the right side of each panel. 


was previously considered ‘natural’. While a great deal of attention had been focused on 
relatively simple models such as the CMSSM or mSUGR.A (with four and three parameters 
respectively) or the NUHM1,2 with five and six parameters, pure gravity mediation models 
can be described with as few as two parameters at the cost of a mass spectrum which 
approaches the PeV scale. As we hope the actual theory of nature is in the realm of 
experimental science, it is imperative to find means to test these models. Here we have 
examined one additional possibility for testing these models despite their generally heavy 
mass spectra. 
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PGM theories, with all their economy, are still able to resolve many of the questions 
their lower energy cousins (such as the CMSSM) were motivated from. These include 
the ability to achieve gauge coupling unification at the GUT scale, radiative breaking of 
the electroweak symmetry, the stability of the Higgs potential, and they also provide a 
suitable candidate for dark matter. The latter is definitely more difficult in PGM models, 
as the wino is usually the lightest supersymmetric particle and as such would require 
a wino mass near 3 TeV to supply the correct relic density. This pushes the gravitino 
mass up to several hundred TeV. Alternatives within PGM are possible if /i is relatively 
small and the Higgsino is the lightest supersymmetric particle [58] or if the theory contains 
additional vector-like states and bino-gluino co-annihilation controls the relic bino density 
[59], or even axion dark matter [58, 60]. In contrast to their lower energy counterparts, 
PGM models have a relatively easy time obtaining a Higgs mass in agreement with the 
experimental measurement [5]. 

Thus experimental verification of PGM models remains challenging. While there is 
the chance that the lightest supersymmetric particle is within reach of the LHC, the bulk 
of the PGM spectrum is not. Here we have calculated the proton lifetime in PGM models. 
We have found that typically the lifetime is long and in many cases significantly above the 
current experimental bounds. However in cases where m 3 / 2 is relatively small and tan f3 is 
relatively high, the proton lifetime is low and may be at the level of current experimental 
searches. While proton decay itself, can not point directly to PGM supersymmetry, it 
may provide one more handle on an ever increasingly elusive theory beyond the standard 
model. 
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Appendix 

A Minimal SU(5) Notation and Conventions 

Here, we review the minimal SUSY SU(5) GUT [29, 30] and clarify our notation and 
conventions. In these models, the MSSM matter fields are embedded into a 5 and 10 
representations of the SU(5) gauge group for each generation. Let and Jh be the chiral 
superfields in the 5 and 10 representations, respectively, with i denoting the generation 
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index. These fields decompose into the MSSM superfields as 
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where a = 1,2,3 denotes the color index. The MSSM Higgs superfields, on the other 
hand, are embedded into a 5 and 5: 
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where the last two components are the MSSM Higgs superfields, 

h r »=( 1 ) • h >= (h!-) ■ < 33 ) 

The other piece of the 5 and 5 Higgs bosons are labeled by H g, and Hca and will be 
referred to as the color-triplet Higgs bosons. 

The gauge boson of SU(5) is a 24. In supersymmetry this corresponds to a real vector 
superfield, V' 4 , where A = 1 ,... ,24 represents the gauge index. V A can be decomposed 
into the SM gauge fields, plus the additional massive gauge bosons of SU(5) breaking, as 
follows 


V = V a T a 
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where T A is the generator of the fundamental representation of the SU(5), and G, B, 
and W denote the MSSM gauge vector superfields with there associated generators. The 
massive gauge bosons associated with the breaking of SU(5) typically referred to as X a 
and Y a will be called just the X-bosons with definition 
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Figure 7: One-loop Higgsino-exchanging diagram which gives rise to the dominant con¬ 
tribution to the baryon-number violating four-Fermi operators. Gray dot indicates the 
dimension-five effective interaction, while black dot represents the Higgsino mass term. 


Here a,/3,... denote the SU(2)i indices. 

The simplest means of breaking SU(5) to the SM gauge symmetries SU(3)c-<g)STJ(2) 
U(l)y is via an adjoint 24 discussed in the text. The 24 decomposes as follows: 


£ = £ A 


rjpiA 




1 ( 2 
+ 2\/l5 \0 



(36) 


Without losing any generality and for simplicity, we assume all SU(5) breaking occurs 
along the £24 direction which is separated in the above equation. 


B Proton Decay 

In this appendix, we give additional details of our calculation of the proton lifetime. The 
important Wilson coefficients arising from integrating out the colored Higgs triplet are 
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These coefficients are then evolved down to the SUSY scale using 
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where / = 1, 2 and Q is the renormalization scale 
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At the SUSY scale Qs, the sfermions are integrated out via the diagram in Fig. 7 to 

Cf = C l e abc (u a R1 d b R J(Ql 3 -L L3 ) , (39) 

ftfr 

i 

where i — 1,2 and 
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These Wilson coefficients C\, which are initially defined at the SUSY scale, are then 
run down from the weak scale using [44] 
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At the weak scale the Lagrangian takes the form 

£(p ->• K + D t ) = C usd [e abc {u a R s b R )(d c L v r )) + C uds [e abc (u R d R )(s c L u T )\ , 

with 

C U sd = ~VtdC2{ r nz ) , 

C uds = —VtsCiirriz) ■ 


(42) 
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(44) 


The new Wilson coefficients C usd)U ds are then further run down to the hadronic scale 
Qhad = 2 GeV. Below the electroweak scale, the RGEs of the Wilson coefficients are given 
by 
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at the two-loop level [61]. The solution for this equation is 
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This long-range renormalization factor is computed to be A R = 1.247 and appears as a 
multiplicative factor to the Wilson coefficients defined at the weak scale. The Wilson 
coefficients at the hadronic scale are then 


C U sd,uds(Q had) C usdud s\m Z ) A R . 

The partial decay width for p —y K + u is then found to be 


T(p ->• K + u ) = 


m p 
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(48) 


where m p and m^ are the proton and kaon masses, respectively, and 


A(p ->■ K + v) = C usd (Q had )(K + \(us) R d L \p) + C uds (Q had )(K + \(ud) R s L \p) . (49) 
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